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ABSTRACT 


The  expressions  for  kinetic  and  strain  energy  for  longitudinal  stress 
waves  in  a bar  are  considered,  first  in  a one -dimensional  model  in  which 
cross  sections  are  assumed  to  remain  plane;  and  stresses  in  the  radial  and 
circumferential  direction  are  assumed  to  be  zero.  From  this,  an  equation 
of  motion  is  derived  which  is  used  to  determine  the  speed  of  longitudinal 
sinusoidal  stress  waves  as  a function  of  wave  length.  Secondly,  a simpli- 
fied three-dimensional  model  is  considered  where  the  axial  motion  is  a 
parabolic  function  of  the  radius,  from  which  the  speed  of  sinusoidal  stress 
waves  is  derived.  The  derived  expressions  are  compared  with  previously 
published  solutions. 


INTRODUCTION 


In  the  simplest  theory  of  longitudinal  stress  waves  in  bars,  in  which 
only  axial  stress  and  strain  is  considered,  the  following  equation  is  ob- 
tained: 

E = pLi  . (1) 


Rayleigh^  considers  a correction  to  take  into  account  the  kinetic  energy  of 
lateral  motion,  which  he  states  as: 

L a 


■1 1.  v^  d X r d r 
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He  then  assumes  that 


V = 


au 

- V ^ r 
3x 


(3) 


This  is  true  if  = <Tq  = 0,  for  then  from  Hooke’s  law^ 


du 


and 


‘ = 
r ar 


IT 


au 

dx 


which,  by  assuming  that  du/dx  is  not  a function  of  r (i,e,  plane  waves)j 
integrates  to; 


du 


V = -vr  ^ + cj  . 


But  when  r = 0,  v = 0 and  therefore  c j = 0,  Thus  we  obtain  expression 
(3),  Differentiating  (3)  with  respect  to  t we  get; 


^ ^2 

av  a u 

— m u r - — 

at  at  ax 


(4) 


If  we  substitute.  (4)  into  (2), 

.L* 


= .4  [ 


T2  = W p \ i r^ 


»2 
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dx  5t , 


r d r d X , 
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We  next  integrate  with  respect  to  r 


T = 
2 


= IT  p 

^O 


- 2 4 / 2 

V a I a u \ 


at  dx 


2 2 2 

For  3 rou  nd  bar,  it  a = A,  and  a =2k  , then; 
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The  correction  factor  for  the  frequency  of  natural  vibrations  leads  to 
the  following  expression  for  the  wave  speed  c of  a longitudinal  sinusoidal 
stress  wave  when  a/^  ^ 1 : 

Love^  used  the  same  energy  correction  for  lateral  inertia  (5)  to  obtain 
a differential  equation  of  motion  for  longitudinal  compression  waves  in  bars: 


P 


2u2 
V k 


T7T?)  ' 


(7) 


From  Love's  equation,  it  may  be  deduced  that  for  longitudinal  sinusoidal 
4 

wave  s. 


c 


E 1 

P 1 I *>  4 i > a vZ 
' 1 + 2n  V 


(8) 


This  is  plotted  as  curve  II  of  Fig.  1.  When  ^ 1,  (8)  reduces  to  (b). 

Thus,  Love's  equation  predicts  that  the  shorter  the  wave  length,  the  lower 

5 

the  speed,  which  is  not  in  accordance  with  observed  data,  nor  with  the 

4 

Pochhammer  exact  solution,  which  is  plotted  as  curve  I of  Fig.  1.  The  exact 
dependence  of  wave  speed  c on  wave  length  s not  in  closed  form,  so  that 
numerical  results  are  not  easy  to  obtain,,  nor  is  the  solution  in  a form  for 
solving  problems  of  impact,  etc. 

It  will  be  shown,  that  with  the  Rayleigh  assumptions,  we  may  add  to 
the  expression  of  strain  energy  a term  representing  the  shear  strain  energy 
of  the  bar  which  will  give  an  improved  estimate  for  the  velocity  of  propaga- 
tion of  longitudinal  sinusoidal  stress  waves. 


THE  ENERGY  EXPRESSION  FOR  THE  ONE-DIMENSIONAL  THEORY 

We  utilize  the  method  that  Love  used  to  derive  (7);  that  is,  obtain 
complete  expressions  for  T and  V that  are  consistent  with  the  assumptions, 
and  vary  (T-V)  to  obtain  the  differential  eqviation  of  axial  motion.  We  make 
the  same  assumptions  that  Rayleigh  made;  that  •'’j.  = «q  = plane 

sections  remain  plane,  or 
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(a)  V = -vr  — , 


u = u(x,  t)  . 


Now, 


V - dv.au 
r 3x  5F  ■ 


But  9u/3r  = 0 from  9(b),  and  using  9(a)  , 


Yxr  = rr 

ox 


The  strain  energy  of  shear  in  a small  volume  dp  is  then 


^^2  = 


xr  ^ G Z G 2 2 / d u I 

T—  = T , dp  = ^ V r y dp  . 


If  we  integrate  over  the  entire  bar. 


^ , !\Z  r-  2 4 

,,  \G22du\,  j ,_Girva 
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The  strain  energy  of  longitudinal  compression  is; 


= I r + ^ j,  (I:) 

v P o 


The  kinetic  energy  of  longitudinal  motion  is: 


L.  a 


J J 


r rl  X = 


^ 7 

r , , 2 


(tr 


.«•  « . J 


-5 


or 


T,  = 


dx 


(13) 


The  variational  equation  of  motion  is 


t 


or 


6 dt(T-V)  = 0 , 
'b 


c 


(14) 


By  substituting  expressions  (5),  (11)»  (12),  and  (13)  into  (14),  we  obtain; 
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dx 
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and  integrating  the  terms  in  (15)  by  parts, 
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If  we  make  the  indicated  substitutions  into  (15)  and  equate  the  coeffi- 
cient. of  8u  to  zero,  we  obtain  the  differential  equation  of  axial  motion: 

0^u  , 2,  2 d^u  . ,,  d^u  , 2 2 d^u  „ 

- p — y + p V k — X y + E — j - G k V — X = 0 • (lb) 
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ANALYSIS  OF  ONE -DIMENSIONAL  SINUSOIDAL  STRESS  WAVES 

To  find  the  relationship  between  the  phase  speed  c on  the  wave  length 
A.  consider  a sinusoidal  wave: 

u = Dcos(2x  - 0#t)  (17) 

and  substitute  (17)  in  (!&)  to  obtain: 


2 , 2,  2 „ 09^  _ . 2 2 

p 0)  + p V k — -j  — E — Y — G k V 

c c‘"  c 
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p + p v“k“  ^ - il-  - Gk^v"  = 0 


Now,  substituting  oi  = Ivc/^  into  (I7a): 


. 2.2  4n^  E „.2  2 4if^ 

p k 


(17b) 


Since 


2(1  + v)  * 


we  obtain 


. . v^k^47T^,  _ E . k^v^2ir^x 

) - -p.  • 


(I7c) 


z z 

Finally,  for  a round  bar,  k = a /2  , therefore 


1 + 2v 


2 2/  \2 
V ir  / a j 

TTT  I ji 
, z z/iV 


Curve  III  of  Fig.  1 is  a plot  of  expression  (18)  for  v = 0,  29. 

With  Eq.  (18),  the  asymptotic  value  of  c as  a/7\.-^oo  is 


c = 


2 2 
V IT 


2v  ir  (1+v) 


E 1 
P 2(1  + v) 


I 

= \|  £ = 0.621  - , 

IP  P 


-■  't- 


which  is  the  wave  speed  of  distortional  waves  in  isotropic  media,  but  it 
should  be  the  speed  of  Rayleigh  surface  waves,^  which,  for  v = 0.29,  is 


= 0.9258  = 0.575  ^ f . 


The  lack  of  fit  of  the  L.ove-Rayleigh  carve  ii  to  the  Pochhammer  curve 
I is  due  to  the  lack  of  consideration  of  shear  strain;  while  the  lack  of  fit  of 
the  derived  curve  III  is  due  to  the  inaccuracy  of  the  two  assumptions  by  which 
(18)  was  obtained.  It  i?  obvious  that  assumption  (9a)  is  inaccurate  because 


j.  --  u.  y ■ 


^8  " 

tha  radial  and  circumferential  stresses  are  not  zero  for  short  wave  lengths; 
and  (9b)  is  inaccurate  because  the  axial  displacements  of  waves  of  short 
wave  length  are  larger  near  the  surface,  but  the  latter  assumption  requires 
plane  sections  to  remain  plane, 

7 

Mindlin  and  Herrmann  have  obtained  an  expression  for  c that  con- 
siders the  effect  of  shear  strain  for  a one-dimensional  theory,  and  this  is 
plotted  as  curve  IV  in  Fig.  1.  They  obtained  their  results  by  expressing 
Newton's  law  in  the  radial  and  longitudinal  direction  in  terms  of  u and  v. 
By  assuming  (i)  v is  proportional  to  r and  (ii)  plane  sections  remain 
plane,  and  then  integrating  the  expressions  over  r , two  coupled  equations 


are  obtained  in  terms  of  u and 


which  arc: 


-2 

,2  ,,2  ^ a,,2..  , , „2,  du  2 •• 

a K p y-  - 8K,  (Krp)-4aK,  X.  — = p a v 


1 


3v  2 - ^ 

2 a X - + a (X  + 2 p)  ^ - p 3 u ; 

8x 


3x 


(20) 


X and  p are  Lame's  constants;  K and  K,  are  abitrary  constants  intro- 
duced to  improve  the  fit  of  curve  IV;  and  Vj=  (a/r)v  , Then  Vj  is  eliminated 
from  the  two  expressions  to  obtain  the  following  equation,  containing  only  u ! 


E 3^u 


3x 
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8K“(x  + p)  ar 


2 ^2 
a K u 

8Kj^(X  + p) 


2 

3x 


X-f 

P 


X+  2 p 3^u 

. P a? 


a^u 


3x 


g 

g 


(21) 


= 0 


The  discrepancy  between  the  Mindlin  and  Herrmann  curve  IV  and 
curve  I can  be  attributed  to  the  assumption  that  plane  sections  remain  plane, 
and  to  an  inaccuracy  in  satisfying  the  boundary  conditions.  If  expression  (17) 
is  substituted  into  the  second  of  Eqs,  (20)  we  find  that: 


V-' 

•»:•■ 

i,' 

»'■ 


f 

b- 
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= -t^  D ( 1 4 P-— •)  sin  ~ (x-ct)  . 


A 


(22) 


Since 


\ / V , 3v  , 3u . , - dv 

r = \(—  + -5-  + •5—)  + 2 UL  ^ j 

r 'r  3r  3x'  ^3r* 


by  substitution  we  obtain 


0-^  = (6|j.-  2pc^  + D sin  (x-ct) 


(23) 


or 


£/E  2.  TT_,.  2tt/  .V 
( — - c ) -T-  D sin  - r-  (x-ct)  . 
r V ' p 'A  A 


(24) 


Thus  is  not  zero  on  the  lateral  sui  faces  of  the  bar;  nor  is  • 


Mindlin  and  Herrmann  state  that  if  the  Rayleigh  assumption  (9a)  is  put 
into  (21),  it  reduces  to  (7).  Equations  (20)  also  include  (16)  if  condition  (9a) 
is  substituted.  From  (9a) 


3u  r 
53E  = a ''1 


so  that 


V j = - V a 


du 

5^ 


(25) 


If  we  substitute  (25)  into  (20)  and  indicate  derivatives  by  dots  and  primes, 

■> 

3 


4 7 «t>  2 r . , T ’ 3 

- a'’  K*  pvu  +4IvjajX,(2v-l)+2vpJu  =vpa  u 


a^(-2X,v  + X+2p)u 


2- 

= p a u 


(26) 


Now 


, ^ _ vE(2v-i)  , 2vE  _ „ 

X.v2v-l)  + 2vp=  + i(T+T)  - ° 


(27) 


and 


2Kv  4\+2p  = X(1 -2v)4  2p  = + ^ (^^^■  v)  = ^ ‘ 


S 


*r  - . 'V 
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SECOND  APPROXIMATION 


u = u (x,  t)  (1  -p 


Then 


„ ./du.dv.v^,,  dv 

0 = M-5T:  + ^ + r + 2^  arr 


Now  substitute  £ur  u from  (31), 


If  we  substitute  (27)  and  (28)  into  (26)  we  get 


,,2  » 

K ^vu  = vpu 


E u = p u 


The  second  of  these  is  Eq.  (1).  If  the  first  of  these  equations  is  multiplied 
by  V k^>  then  differentiated  with  respect  to  x,  and  added  to  the  second 
equation;  we  obtain 

-pu+pv^k^u  +Eu  - p.  k^  V u = 0 . (30) 

Now,  if  K is  set  equal  to  unity,  (30)  reduces  to  Eq.  (16), 


The  assumption  that  cross  sections  remain  plane  as  longitudinal  stress 
waves  pass  is  physically  unrealistic,  lor  Eq.  (10)  gives  a finite  value  of 
shear  strain  and  stress  at  the  surface  of  the  bar.  To  eliminate  these  shear 
strains,  it  is  necessary  that  the  axial  displacement  be  a function  of  the 
radius  as  well  as  of  x and  t.  L.et 


and  assume  to  be  zero.  This  determines  v since 


r = X(«  + t + €-)  + 2u  € 
r ' X r 8'  r 
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Solving  (34)  for  v- 


o / p r I 


Since  v=Q  when  r = 0,  Aj=0.  To  deterniine  p,  we  use  the  condition 

that  the  shear  strain  at  the  surface  is  zero: 


“ 9r 


r 2 ® 

a 9x 


“o  / Pr’  1 


T « -20  v®“o,  vPa 

Y I = 0 = •~*T  XI  *f  ; jA'*'  J 

Jr=a  a o a 


Hence 


2 ^ '^o 


2u  - 


r 


For  sinusoidal  stress  waves  • (17),  the  value  of  p is: 


2.  4v  /ira\ 


'i';:  I 


The  strain  energy  of  compression  is 


^1  " i J ^"‘x*x+'"r*r  + ”’9 


s*<r 


where 
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We  next  integrate  over  r and  substitute  (31)  to  get 

V.  f"'sin2fe\x-ct)dx  . 


\A/  ^ 


\J^I 


whe  re 


♦ j = 1 - P + 0.  3439  P . 

The  strain  energy  of  shear  is 
Li  a 


^2  = 


? 1 1 


Y Zitrdrdx 


o ''o 


Since 


L 


I 

I 

I 


.1 


, = -2fli+  It  . 

” [ 7 \ A / \ 


(3-2v).'=y 


D cos  1^;^  (x  - C t) 

A 


1^2  ^ 

, = GAv^l  — 1 ^ r cos^  ^ (x  - c t)  d X , 

" \M  ^ ^ J„  A 


V, 

i 

whe  re 

A = 1 _ 4 y P . + P 

♦Z  ‘ T (3  - 2 v) 


I • (a, 


2(3  » 2 v)  V w 

+ P" 

4 V IT  \ / 


(‘  • I ttAt}] 


The  kinetic  energy  of  longitudinal  motion,  using  (13),  is 


Integrating  with  respect  to  r. 


sm^  — (x-ct)dx  , 
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' . r*Y 

i' 1/ 

•H 


I 

r 

I 

(48) 

I 

I 

I 

(49) 

(50) 

I 

, (51) 

(52) 

(53) 

i 

i 

I 
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(54) 

(55) 
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where 


= 1 - P + 


The  kinetic  energy  cf  lateral  tnotion,  using  (2)  is 


’■■'I  {■ 


2 2 « r’  . 

- p -j.  + 


(3-2  v)^  ? \ 


^ 2TT  r 


d r d X 


If  we  integrate  and  substitute  (31) 


T,  = pv‘ 


where 


A L 

2 . ^ „2  r 2 2tt  , . 

I c A ^2^  J vx-c  t;  dx 


»*  - 1 ^ K I 

'-2  ■ ^ ■ T TT-^ 


2(3  - 2v)“ 


We  perform  the  variation  on  T - V;  and  obtain 

p(^)  ASjDtpv^^^l  A ^ {jD-Ea/^I  ■tjD 


from  which 


e _ (jif 

7p  S'  4.  7 

' (ai 


which  has  the  same  form  as  Eq.  (18).  Expression  (6l)  has  been  plotted  as 
Curve  V in  Fig.  1,  where  it  is  seen  that  it  fits  the  exact  curve  I most  closely. 
The  discrepancy  between  I and  V is  due  to  the  inaccuracies  of  assuming  a 
parabolic  variation  of  u with  r,  and  assuming  = 0.  Figure  2 is  a 
comparison  of  expression  (31),  the  axial  displacement,  with  the  exact  solu- 
tion  as  found  by  Davies  for  aA  = 0.  196.  The  radial  stress  of  the  exact 
is  hIso 


"r 

•'  ^ - .'T 
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CONCLUSION 

The  derived  one -dimensional  theory  of  axial  motion,  expression  (16)^ 
may  be  useful  in  a qualitative  way  say  for  determining  the  dispersion  of 
stress  waves;  but  since  this  equation  is  based  on  the  assumption  that  there 
is  no  radial  stress  and  that  cross  sections  remain  plane,  for  greater  ac- 
curacy the  Induced  curvature  of  the  cross  sections  should  be  taken  into 
consideration. 
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SUMMARY  OF  SYMBOLS 


a 

A 

c 

D 

E 

Q 

G 

k 


- radius  of  bar 

- cross-sectional  area  of  bar 

- phase  velocity  of  longitudinal  stress  waves 

- amplitude  of  axial  displacement 

- Young’s  modulus 

- elongation 

- shear  rigidity,  = E/2(l+v) 

- polar  radius  of  gyration 


K,  - arbitrary  constants  in  Mindlin  and  Herrmann  equation 


L. 

r 

t 

T 

T. 


u 

V 

V, 


- length  of  bar 

- radial  coordinate 

- time 

- kinetic  energy 

- kinetic  energy  of  longitudinal  motion 

- kinetic  energy  of  lateral  motion 

- displacement  in  x (longitudinal)  direction 

- displacement  in  r (radial)  direction 

- strain  energy  of  longitudinal  strain 

- strain  energy  of  shear 

- longitudinal  coordinate 


Y 
A 
0 
A 

V 


- shear  strain 

- dilatation  = € +«_+*, 


X r 0 

- coordinate  in  circumferential  direction 


- wave  length  of  sinusoidal  waves 

- Lame’s  constant,  = vE/(1+v)(1  - 2v) 

- Lame's  constant,  = G = E/Z(1  +v) 

- Poisson’s  ratio 

^2“  correction  factors  in  expression  (61) 
p - mass  density 

(T  - normal  stress 

•v  - shear  stress 

^ - volume 

correction  factors  In  expression  (61) 

A Lt 

0)  - 2n  times  the  frequency 
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